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A model study of the singlet excitons in the Ceo molecule 
with emphasis on the Coulomb interaction between excited 
electron and hole leads to a physical understanding of the 
interaction effects on the absorption spectra and to a new 
identification of the forbidden excitons in the third-harmonic 
generation spectra. These conclusions may be tested experi- 
mentally on the model predictions related to the optical Kerr 
effect. The model shows that, with sufficiently strong inter- 
atomic than onsite interaction, a T2G exciton could have very 
low energy or become unstable against the closed-shell ground 
state. Properties of these interesting cases beyond the Ceo are 
briefly examined. 
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I. INTRODUCTION 

Since the discovery of Cqq%, several experimental and 
theoretical studies have suggested that undoped and 
doped solid 0*60 are strongly correlated electron systems. 
The observational basis includes the unusually high su- 
perconducting transition temperature of the alkali-metc 
doped fulleritescl, the existence of soft ferromagnetisr 
and the strong Coulomb interaction effects in the Auger, 
direct and inverse photoemission spectraa. The theoret- 
ical motivation is based on the narrow bandwidth com- 
pared to the strong intra-molecular interaction. This nat- 
urally leads to theories based on the Hubbard model with 
onsite interaction on every carbon atomoEl. Strong corre- 
lation is used to find a mechanism for superconductivity!! 
Parity doublets of the lowest unoccupied molecular or- 
bitals (LUMO's) of the Cqq molecule are used for 
the pairing mechanismcl. Coulomb interaction of the 
electrons within a single Cqq molecule has also taken 
into account through configuration interaction by sev- 
eral quantum-chemistry calcuk^ionsBilM On the other 
hand, quasipaptiek: correctionlljo to the local density 
approximatioriialia leads to the conclusion that, while the 
correlation effect is large, the electronic structure, in Cqq 
is nonetheless that of a standard band insulatoiO. Simi- 
larly, the superconducting transition temperature of the 
alkali-doped-jftio has been explained by the usual phonon 
mechanismll2Hl3. 

A common feature among the theories mentioned 
above is the assumption of the closed-shell ground state 
for Cqq with the set of highest occupied molecular or- 



bitals (HOMO's) of symmetry h u and h g and the set of 
LUMO's of symmetry t% u and t\ g . We ask the question 
whether such a ground state is stable against the excita- 
tion of an electron from a HOMO to a LUMO. If not, the 
new ground state would lead to low-lying excited states 
of a quite different nature. In this paper, we use the Cqq 
molecule as a paradigm for molecular solids to investi- 
gate the stability of the closed shell ground state. By 
using a simple model, we hope to understand the factors 
governing such instability and to explore consequences, 
such as in nonlinear optical properties. Even if the ex- 
cited states are only low in energy without causing any 
instability, they could play an interesting role in some 
properties, especially superconductivity. The method of 
this study may be applied to molecular solids and quan- 
tum dots besides Cqq. 

First, we address the energy ordering of the closed- 
shell state and, the one electron-hole pair excited states. 
Shirley et al.E3 used a molecular orbital model to make 
a comprehensive study of the exciton energy spectrum 
in solid Cqq. We adopt a similar approach but use fur- 
ther simplifications. We restrict our attention to a single 
Cqq molecule since ixuthe solid the weak overlapping be- 
tween the moleculeso would not qualitatively affect our 
results. We use a nearest- neighbor Jtight-binding model 
for the one-electron it otbitalsElO^Ed with the help of 
symmetry considerations^. The essential results of this 
model are given in Sec. |l. In Sec. [II, the energies of 
the electron-hole pair states relative to the closed-shell 
state are determined in terms of the single-particle ener- 
gies and the Coulomb interaction. The interaction terms 
include direct electron-hole attraction and the exchange 
counterparto, interaction on the same carbon site as well 
as between any pairs of carbon sites in the same Cqq 
molecule. The dependence of the pair-state energies on 
three model parameters, the nearest-neighbor hopping 
energy V, the intrasite interaction U (including both the 
direct and exchange contributions) and the typical long- 
range interaction term e 2 /(ei?o) taking into account the 
intersite screening effect e and setting the distance scale 
at the radius of the buckyball, Rq (w 3.5 A) is discussed. 
We find that a scenario of either the closed-shell state 
or an electron-hole pair state being the lowest energy 
state is possible for reasonable values of the three pa- 
rameters. The consequences of both scenarios on the lin- 
ear and nonlinear optical properties are then considered 
in Sec. Our calculations of the linear optical spec- 
tra are compared with experiment and our calculations 
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of the nonlinear properties are used to suggest measure- 
ments which will clarify the situation. Sec. [v| discusses 
the results of our work with regard to the nature of the 
ground state of Cqq. Our tentative conclusion is that the 
comparison of the theoretical and experimental linear op- 
tical spectra favors the closed-shell state as the ground 
state. However, there are some unsatisfactory features. 
We believe that the electro-optic and nonlinear optical 
measurements can clarify the situation. We speculate on 
the possibility of molecules with an excitonic state as the 
ground state. 



II. SINGLE-PARTICLE EXCITATIONS 

Since one objective of this work is the study of the 
stability of the closed-shell ground state against the low- 
lying excited electron-hole pair states, the first order of 
business is to construct the most relevant one-electron or- 
bitals in C 60 , namely the HOMO's and LUMO's. Their 
approximation by 7r orbitals seems to be well established 
for low-energy excitations. For instance, the weights of 
the radial orbital for the h u and t\ u states were about 
98% and 95%, respectively, according to Laouini et al.El 
In the discussion of the electron-hole interaction the small 
cr-contributions may be negligible. The tight-binding 
Hamiltonian with nearest-neighbor hopping may be built 
from the molecular ir orbitalaH3: 



km P (x) = N l^2 e i^p{s)X ( x - #og) 



(2.1) 



where x( x — RoS) is the component of the p wave func- 
tion centered around the atomic sites i?og pointing along 
the radial direction. We have neglected the difference in 
bond lengths of the inequivalent bonds. Ni gives the nor- 
malization of the molecule state. The quantum-jiumber 
I is used to index the irreducible representationsEj of the 
icosahedral group Ih , the symmetry group of the bucky- 
ball, namely a, t\, h, t2, and g, with degeneracies 1, 2, 
5, 3, and 4, respectively. The quantum number m runs 
over the degenerate states of each irreducible representa- 
tion. The quantum number p denotes the parity of the 
state. Its introduction indicates that the full symmetry 
group of 6*60 is Ih x Z2, where Z% is the two-element 
group consisting of the inversion operator and the iden- 
tity. The coefficients for different sites may be related to 
each other bya 
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standard transformation matrices in a rigid body, the so- 
called Wigner-Z? functionsEH. 

The irreducible representations D l mm , (uj g ) of the coef- 
ficients e; mp (g) reduce the nearest- neighbor Hamiltonian 
in units of the happing parameter — V to a set of Hamil- 
tonians given byu: 



Yd 1 ' , 

/ j mm 



(2.3) 
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^(sin(26 ),0,2 + cos(20 o )), where the angle 
cos _1 (l/v / 5) is defined by the geodesic arc be- 



where the i-sum only runs over the three nearest neigh- 
bors fi of the site e. The spin degrees of freedom are 
understood. For readers interested in generating the ir- 
reducible representations^], we record the coordinates 
e = 

2©o 

tween two neighboring vertices of the icosahedron. The 
rotations to the nearest neighbors, fi, from the atom at e 
can be given by the Euler angles (a — 0, f3 — 20o, 7 = tt), 
(a = 2tt/5,/3 = 0,7 = 0), and (a = -2n/5,/3 = 0,7 = 
0), respectively for i = 1, 2, 3. 

The LUMO's and HOMO's of interest belong to the 
representations t\ and h respectively, which are isomor- 
phous to the spherical harmonics I — 1 (p wave), and 
I = 2 (d wave) . The state degeneracy is then 21 + 1 and 
the normalization Ni = \J (21 + l)/60. The correspond- 
ing Hamiltonians for these states are 
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where cos(29 ) = l/\/5 and cos(^) = |(-1 + y/E) are 
used. 

The eigenvalues A of the reduced Hamiltonians for the 
two LUMO and two HOMO levels are, for I = 1. 



A i+ = i(-3 + V5), 
Ai_ = i[(3 + V5)/2 



\J(19-V5)/2], 



(2.5) 



where oj g is the rotation bringing the radial vector from and for 1 = 2, 
atomic site e (see Fig. 1, not to be confused with the 
coefficients e; mp (g)) to site g. For the regular three- and 
five-dimensional representations (/ = 1, 2) under consid- 
eration the (21 + 1) x (21 + 1) matrices D l are simply the 



A 3 _ = i(-1 + V5), 



(2.6) 
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in agreement with other calculations0€3O. These or- 
bitals are associated, respectively, with the symmetry 
hg, hu, h u , and h g . The single-particle energies of these 
molecule states are 



-Xi P V, 



(2.7) 



where Ip runs over the indices 1+ (for the representation 
tig), 1- {h u ), 2- (h u ), and 2+ (h g ). 
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For these parity double ts w ith not so very different en- 
ergies, the Hamiltonians (2^4) give the normalized eigen- 
vectors with (21 + 1) components (with x = Ai_ — 2 

and tany = \/bx/2) at carbon site e as 

1 / 1 \ i / - sm ^ 

ei+(e) = -j= I I , ei_(e) = -j= ( x/2cosy 
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triplets for the LUMO states and quintets for the HOMO 
states. The complete eigenvectors with the components 
for the other atoms follow by rotation (2.2). Usiftg this 
equation, the definition of the Wign er-D functionsEH, and 
the form of the vectors in Eq. ( |2.8| ) one can easily show 
the parity of the states, e im ±(-g) = ±e Zm ±(g). The 
pairs of vectors for the atoms at the sites e and — e give 
instructive examples. The corresponding transformation 
matrices D l mm ,(0,ir,Q) = (-l) l+m 8 m - m ' gives eigenvec- 
tors at — e, which fulfill the parity condition. Since the 
rigid-body transformation from e to — g may be related 
to a product of transformations from e to g and e to 
— e, the above property is also valid for arbitrary atomic 
positions g. 

Fig. 2 shows schematically the two LUMO levels and 
two HOMO levels and their associated states. The e nerg y 
scale is set by the hopping matrix element V in Eq. (2.7). 
The single-particle energy difference (ei_ — £2-) is taken 
to be 3.5 eV between hu and h u peaks in solid Ceo mea- 
sured by the, phqtoemission and inverse-photoemission 
experimentsa'EScj. This yields an estimate of V = 4.626 
eV. The brackets V = (3.83, 6.61) eV represent the uncer- 
tainty of this estimate. The lower value arising out of the 
finite band widths is taken to be the midpoint between 
the band onset at 2.3 eVllI and the peak-to-peak difference 
at 3.5 eV. The higher value of V comes from the estimate 
of 5 eV as the difference between the electron affinity level 
and the ionization potential of the Cqq moleculecJ. The 
three values of V yield the single-particle excitation spec- 
trum £i+ = 4.63 (3.83, 6.61) eV (t lg ), £i_ = 3.50 (2.90, 
5.00) eV (t lu ), £ 2 - = eV (h u ), and £ 2+ = -1.77 (-1.46, 
-2.52) eV (h g ) with respect to the position of the highest 
occupied state. The discrepancy between these values of 



V and the LDA derived V = 2.72 eV0 represents the 
phenomenological fit of the former to the renormalizcd 
one-particle energies so that, when the interaction be- 
tween two single-particle excitations is considered later, 
the one-particle energies should not be further modified 
by the interaction. 



III. EXCITONS 
A. Symmetry-adapted electron-hole pair states 

In this paper, we shall consider only electron-hole ex- 
citations without spin flip, i.e. only singlet excitons. We 
have calculated the energies of the triplet excitons, which, 
devoid of exchaiige interaction terms, lie slightly lower 
than the singletslij. The possibility of magnetism involv- 
ing the triplet state will be left to a future study. 

Consideration of the electron-hole pair excitations of 
the 7r-electron system of Cqq, which are lowest in energy, 
can be restricted to the level scheme of Fig. 2 with the 
empty levels hg and hu and the occupied states h u and 
h g . The pair excitations contain products of the type 
hp e x h Ph with the single-particle parities p e ,Ph = +1 
(g) or = —1 (u). With the pair parity P = p e -p^, where 
P runs over the same values +1 (G) and pi 1 (U) as p e 
and ph , the pair states have the symmetry^ 



hp e x h Ph 



Tu 



T, 



2P 



(3.1) 



That is, each of these four 15-dimcnsional product 
representations of a singlet electron-hole pair from 



ti Pc splits up into 



LUMO/HOMO of the type h Ph 
two three-dimensional representations, Tip and T2P, 
one four-dimensional Gp representation and. one five- 
dimensional Hp representation (P = G, C/)EJ. Among 
them is the dipole-allowed pair excitations T\u for elec- 
tron and hole with opposite parity. For these optically 
observable excitons, there is no need to consider the four- 
fold degenerated a™ hole level which is either somewhat 
below the h g levelE] or degeperate with it withiruthe ap- 
proximations consideredB'r 2 rr 3 l, from the relationE2l 



h Pe x g g — T 2 p + Gp + Hp, 



(3.2) 



which contains no representation T\u of the electric- 
dipole-allowed excitons. Symmetry reasons also dictate 
that there is no configurational interaction between the 
Aq ground state and the low-lying pair states considered. 

We develop a method of computing the excitonic 
states in terms of the symmetry-adapted electron-hole 
pair states (and incidentally gained some physical in- 
sight into these pair states) by exploiting the close re- 
lation of the representations of the symmetry group Ih 
of the Cqq molecule with ihe transformation properties 
of the spherical harmonicsE3 which differ by a small per- 
turbation. The (£i Pc ) of the LUMO and the (h Ph ) of the 
HOMO correspond to single-particle angular momentum 
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states with the quantum numbers I and m (I = 1,2 and 
— I < m < I). The angular momentum addition rules 
would yield the symmetry of the resulting electron-hole 
pair states to be those of the spherical harmonics L and 
M with L = 1,2,3 and -L < M < L. Indeed, the pair 
states have a three-dimensional representation Tip corre- 
sponding to L = 1 and a five-dimensional representation 
Hp corresponding to L = 2. However, since 6*60 does not 
have complete spherical symmetry, the L — 3 states split 
into two groups, a three-dimensional representation Tip 
and a four-dimensional one Gp. The symmetry-adapted 
electron-hole pair spin-singlet states may be written as 
linear combinations 



\LNp e p h ) 



E E 

m e = — 1 mh — — 2 



(3.3) 



where the operator cf pm (ci pm ) creates (annihilates) an 
electron in a molecule state \lpm > with a single-particle 
parity p. The coefficients on the right of Eq. (SjLare 
related to the Clebsch-Gordan coefficients 
(L = 1, 2 with -L < N < L and L = 3 with iV 



,TBby 
0,±1) 



^m e m h (LN) — ( 1) mh C^ 2 _ mh 



C 2 (3,±3) = (-l)-™'* 



c lm e 2-m h 



C£. mfc (3, ±2) = (-!)-"» 



2-m h 



-f3±3 



(3.4) 



The symmetry- adapted pair states are chosen such that 
quantum numbers L and N, with L = 1 corresponds to 
the basis of the irreducible representation Tip, L = 2 
to £T P , L = 3 and N = 0, ±3 to T 2 p and X = 3 and 
N = ±1,±2 to G P . 

The true singlet exciton states are li nea r combinations 
of the symmetry-adapted pair states (3.3) 



\LNPA) = c Ap (LNP)\LNP-p p), 



(3.5) 



where the summation runs over the hole parity p. The 
summation in Eq. (3.5) indicates that pair states of dif- 
ferent single-particle parities may be coupled provided 
that the total parity P is conserved. The fourth quan- 
tum number A labels the two coupled pair states of the 
same symmetry. The four quantum numbers, L, N, P, 
and A span the 60 pair state s (w ithout spin) under con- 
sideration. The eigenstates ([5Jj) of the Frenkel excitons 
are orthonormalized with J2 P c Ap c A'p — ^AA' , following 



the orthonormalization property of the Clebsch-Gordan 
coefficients. 

The symmetry- adapted basis pair states \LNP-p p > 
of Eq. ( p.3j ) with the third quantum number set to p e — 
P-p block-diagonalize the two-body Hamiltonian of the 
7r-electron system including the full Coulomb interaction 
v into 2x2 matrices diagonal in the quantum numbers 
L, N, and P: 



(LNP-p p\H\L'N'P'-p' p') 
^ll'Snn'Spp' {Sp P ' [eip.p — £2 P ] 
+ (LNP-p p\v\LNP-p' p 1 )} , 



(3.6) 



with the electron in the level with the excitation energy 
£i Pc and the hole in the level with S2 Ph and with the 
Coulomb interaction connecting pair states with (p e ,Ph) 
and (— p e ,— Ph)- In our notation system, the exciton 
energy eigenvalues Xltvpa are independent of N for 
L = 1,2, i.e. (2L + l)-fold degenerate. For L = 3, the 
exciton energies for the two symmetry sets of N = 0, ±3 
and of N = ±1,±2 are different but degenerate within 
each set. 



B. Coulomb interaction 



The Coulomb interaction term in Eq. (3.6) includes the 
electron-hole attraction--and the exchange term to the 
electron-hole attractioncil with the diagrammatic repre- 
sentation in Fig. 3. The exchange terms of Fig. 3b only 
exist for the spin-singlet exciton. Assuming non-overlap 
of the p z -orbitals from different carbon sites of the tt- 
like mole cule states, we express the Coulomb term in 
Eq. ( |3~6|) in terms of the single-particle eigenstates of 
Eqs. (EE§) and (O) as 



(LNp e p h \v\LNp' e p' h ) 
15 



J2 Gl N (12p ePh \gg')v(g - g')G LN (12p' eP ' h \gg') 

9,9' 

2 E Gl N (l2p ePh \gg)v( g - g')G LN {l2p' e p' h \g'g'\ 



9,9' 



(3.7) 



with 



G LN (l e lhPePh\SS') 



+lh 

z. E 

e=— h m h =-l h 



Ct[ h mh {LN)e* lM (g)ei hPhmh (g'), (3-8) 



where the coefficients C l ^ m (LN) are defined in Eq. 



4|). The Coulomb potential takes the form 

e 2 1 

v(g) = U5 g0 + — T1 (l-6„ 



J g0), 



(3.9) 
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where U denotes the one-site Coulomb matrix element of 
the p z orbitals. i?o is the distance from a carbon atom 
to the center of the cage and e denotes a dielectric con- 
stant representing the screening of the interatomic (but 
intramolcculc) Coulomb interaction. 



The first term on the right of Eq. (3/7) comes from 
the direct electron-hole attraction and the second is the 
exchange counter part. The factor of 2 may be viewed 
as originating from the spin degeneracy or the structure 
of the singlet. The double sums over the carbon sites 



in expression (3.7) may be reduced to single sums, by 
means of the product relation for the Wigner D-functions 
representing the two rigid-body transformations to the 
sites. Thus, with a simplifying definition for the Coulomb 
matrix element: 

V pp > (LNP) = - (LNP-p p\v\LNP-p' p') (3.10) 

= -^L Fpp ,(LNP) + [H PP ,{LNP) - 2X PP ,{LNP)\ , 

which the N— dependence serves only to differentiate be- 
tween Tip and Gp in the L = 3 case. For L = 1,2 
(-£ < N < L) and L = 3 (N = 0, ±1), the intraatomic 
Coulomb interaction is given by 



F PP ,{LNP) 



(3.11) 
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2L 



T E <?2 M (12P-pp|ee)G iM (12Py p'\ee) 



M=-L 

and the interatomic Hartree and exchange contributions 



M=-L g 



G* LM (12P-p p|ge) 



-e 



(3.12) 



-G LM (l2P-p' p'\ S e), 



X P ALNP)= EE' 



M=-L g 



G* LM (12P-p p|gg) 



— e 



(3.13) 



■G LM (l2P-p' p'\ee), 



where the functions G* la/ a re defined in Eq. (3.8). The 



intraatomic term, Eq. (3.11), includes both the Hartree 
and exchange contributions. 

In Table I are list ed the val ues of these three terms 
evaluated from Eqs. ( 3.11 - 3.13 ). The corresponding re- 
sults using the continuum approximation of Ref. |^ differ 
little for the Hartree contributions, but up to 50 % for the 
exchange terms. Table I indicates that the interatomic 
electron-hole exchange may be neglected in comparison 
with the interatomic electron-hole attraction. This is 
in complete contrast to the intraatomic case, which is 
exchange-dominated since the singlet exciton has twice 
the number of exchange terms of equal magnitude as the 
direct attraction (cf. Fig. 3). Since the prefactor U/60 
is smaller than e 2 /(ei? ), it is evident that the diagonal 



elements, V PP (LNP), in Eq. (3.10) are dominated by the 
interatomic Hartree matrix elements, and are, therefore, 
positive. The off-diagonal elements p p' are strongly 
influenced by the intraatomic exchange. They are often 
negative (except for T\q and Hu). Table I also shows 
that a contact-potential approximation, where the inter- 
atomic Coulomb interactions are neglected, is invalid. 



C. Pair excitation energies and ground-state stability 

With the Coulomb interaction given in Table I, the 
2x2 eigenvalue problems for the 60 Frenkel excitons orig- 
inated from the closest 7r-electron-related HOMO and 
LUMO single-particle states can be solved. For a given 
representation Ti, H, T 2 , or G and a given total parity 
P = ±1, the 2x2 Hamiltonian for hole states of parities 
p and p' may be written as 



H pp > = (E+pA)5 pp , - U(l - S pp >), 



(3.14) 



where 



E = -{E + +E_), 

A = ±(£ + -£_), 

U = V+-(LNP), 
E p = s 1P .p - e 2p - V PP (LNP). 



(3.15) 



The two exciton eigenvalues -Elatpa with A = ±1 and 
the corresponding eigenvectors follow as 



Elnpa 
c AP (LNP) 

where 



e + aVa 2 + u 2 , 

Sa- p> + cosry + S A . p -p ■ sin 77, 



sin(2ry) = U/\/A 2 +U 2 . 



(3.16) 



(3.17) 



The coupling of the Frenkel excitons with the same rep- 
resentation and parity P destroys the simple picture that 
the reduction of the d iffere nce of the single-particle ener- 
gies Eip. p — S2 P in Eq. ( 3.15 ) by V pp defines the binding en- 



ergy of the exciton. The ratio V j\J A 2 + U 2 determines 
the strength of the redistribution of the two coupled ex- 
citons with A = ±1. Therefore, the sign of U plays an 
important role for the actual oscillator strength for the 
excitations of electr on-ho le pairs with different A as will 
be discussed in Sec. 



IV B 



In Figs. 4 and 5 a selected set of pair excitation en- 
ergies are plotted versus the strength of the interatomic 
Coulomb interaction e 2 /(ei?o) for two different values of 
the intrasite Coulomb matrix element U. To avoid clut- 
ter in Fig. 4, only plotted are excitons of representations 
Tip (L = 1) and Hp (L = 2), i.e., a total of eight exciton 
energies with P = ±1 and A = ±1. The other eight exci- 
ton energies for the representations T 2 p and Gp (L = 3) 
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are not plotted since they do not usually appear in the 
optical spectra considered below. Figure 5 compares the 
four lowest pair excitations Tig, T 2 g, Gg, and Hg with 
the Ag Hartree-Fock ground state with even parity. The 
two intraatomic interaction values chosen are U = and 
U = AV. The Coulomb energy U when two electrons are 
in the same atomic p-orbital is usually estimated to be 
U ~ 10-20 eV. If one uses a hydrogenlike wave function 
with an effective nuclear charge z e ff = 3.25, a value of 
U = 17.3 eV resultsa. With reasonable values for the 
hopping parameter V (cf. Sec. |n|), U = AV is close to 
the estimated values. Both Figs. 4 and 5 show that the 
explicitly chosen U value has a minor influence on the 
excitation energies. A strong intraatomic Coulomb inter- 
action U — AV shifts the pair energies slightly to higher 
energies by about 0.05V. 

The dependence on the variation of the interatomic 
Coulomb interaction e 2 /(ei?o) is explored because of the 
uncertainty of the value for the dielectric constant used 
to screen the interatomic interaction. The range goes 
from the limit of metallic screening (e — oo) to the limit 
of the unscreened Coulomb interaction, (e = 1). With 
Ro w 3.5A and the range of the hopping parameter V 
Sec. 0, the maximum value of e 2 /(eRoV) is about 
1.1. In the literature, screening values between 3 and 10 
have been reported. For instance, dielectric constants for 
solid, faillcrites have been determined as e = 3.5, 3.9, or 
4.4c3~E3. A dielectric constant in a model cluster of 7.13 
to 9.86 has been used in Ref. [55| to study the van-der- 
Waals cohesion energy. Other authorso use e = 4.4 and 
6.5_to explain the screening in Cqq clusters. Hansen et 
al.B3 reported a value e = 4.6 derived from a Kramers- 
Kronig analysis of their visible- UV EELS spectrum. For 
e = 4.6, the interaction parameters of e 2 /(eRoV) « 0.19 
(0.23, 0.14) result from V = 4.626 (3.83, 6.61) eV. From 
Fig. 4, evidently the interatomic Coulomb interaction has 
much more influence on the pair-excitation energies than 
the onsite Coulomb interaction. Moreover, its presence 
gives rise to an effective attractive interaction between 
electrons and holes. The exciton binding noticeably re- 
duces the pair excitation energies. Typical reductions 
amount to about 0.8e 2 /ei?o f° r Tig, Hg, Hu, and Gg 
and vanishing intraatomic interaction U . For T\jj (T2g), 
slightly smaller (larger) values of 0.5 (1.0) e 2 /ei?o are 
shown in Figs. 4 and 5. The reductions corresponding to 
finite U values are smaller. 

The drastic reduction of the electron-hole pair excita- 
tion energies by the interatomic Coulomb interaction is 
illustrated in Fig. 5 for the lowest excitations with even 
parity. The figure indicates that for sufficiently weak 
screening, pair energies may even become negative. For 
Tig this happens if e < 6.5 eV/V (when U — 0) or e < 5.5 
eV/V (when U = AV). For V — 3.83 this corresponds 
to e < 1.7 or e < 1.4. For a very strong (i.e. almost 
unscreened) interatomic Coulomb interaction, the lowest 
Tig pair excitation replaces the closed-shell Hartree-Fock 
state Aq as the ground state of the Cqo molecule. This 
transition should be accompanied by a spontaneous sym- 



metry breaking to lift the degeneracy of the lower Tig 
states with respect to the quantum number N. 

The findings in Figs. 4 and 5 are consistent with the 
results of quantum-chemical calculationsEfllj. The en- 
ergies of the 1 Tig, 1 T2G, and 1 Gg multiplets are quasi- 
degenerate, within a range of 0.1 eV, at about 0.5 eV and 
the 1 Hq multiplet is higher, separated by about 0.4 eV. 
Figure 5 (b) (U = AV) indicates a similar situation for 
reasonable parameters of e 2 /(eRoV) w 0.2-0.5 discussed 
above. 

The comparison with the quantum-chemical results 
may be used to fix our model parameters. In the un- 
screened case e = 1 the disagreement is too large. When 
the effect of the interatomic Coulomb interaction is some- 
what screened, e.g. in the case of U = 17.3 eV and an 
intermediate effective screening of e = 4.6, the agreement 
is improved. For the three hopping parameters V = 4.626 
(3.83, 6.61) eV we find 2.86 (2.53, 3.90) eV for T 1G , 3.06 
(2.65, 4.24) eV for T 2G , and 2.75 (2.46, 3.71) eV for G G , 
and 3.16 (2.71, 4.41) eV for Hq- Consequently, we con- 
clude that the effective interatomic Coulomb interaction 
has to be effectively screened. For the hopping parameter 

V = 3.83 eV, e = 4.6 is the appropriate value. For the 
larger hopping parameter V — 4.626 eV the screening has 
to be slightly increased. The largest hopping parameter 

V = 6.61 eV should be excluded since the dielectric con- 
stant becomes unreasonably large. The resulting picture 
is also more or less [Consistent with quasiparticle calcu- 
lations for solid C6oc-3 which yield the lowest excitation 
energy of 2.15 eV for the h u — > t\ u transition. When 
the electron-hole interaction is included!^ this value is 
reduced to 1.57 eV giving an exciton binding energy of 
0.58 eV. This ab initio band calculation also gives peak 
positions for T 2 g, Tig, Gg, and Hg that are redshifted 
in comparison to the peak positions 1.86, 1.94, 2.03, and 
2.30 eV in the fine structure of the forbidden absorption 
edge of the fulleritetS. An inconsistency among the first- 
principles calculations is the energy ordering of the Tig 
and T 2 g excitations. In Refs. [l(],[ll| the symmetry of-the. 
lowest excitation is Tig whereas other calculationsa'Ej 
indicate the T 2 g level to be the lowest one. Figure 5 
shows that the answer depends on the relative strengths 
of the intra- and interatomic Coulomb interactions. For 
[7 = 0, T2G represents the lowest excited state. For 
U = AV, this holds only for e 2 /(eR V) > 0.7. In the 
more interesting region of lower interatomic values Tig 
is favored. For even smaller values of interatomic inter- 
action, Gg becomes the ground state. 



IV. OPTICAL SPECTRA 

A. Optical transitions 

The coupling of light to the 7r-electron system of the 
Cqq molecule is governed by the polarization operator 
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Pn 



ex a v ) c, c v > 



(4.1) 



where x Q is the a-th Cartesian component of the posi- 
tion operator. With the restriction to the LUMO and 
HOMO states described in Sec. [n] and a strong localiza- 
tion of the p z orbitals as in the c ase of the description of 
the Coulomb interaction in Eq. ( |3.9| ), the dipole matrix 
elements take the form (y = Ipm): 



LNp e p h \P a \L'N'p' e p' h 

+ 1 +2 

EV^ / -i \ —mh—mhi r^LN (iV N' 

/ j \ L ) u lra e 2-m h L '1^2-m^ 

m e ,m' = — 1 m h ,mi=— 2 



X{-^PeP'Jm^m' e (^p'h m 'h \ ex a\ 2p h m h ) 

+Sp hP ' h S mh m h , {lpeme \ex a \ lp' e m' e )}. 



(4.8) 



(Ipm \ex a \ I'p'm!) = 

V(2Z + 1)(2Z' + 1) 
eRo 60 



We record here a special case needed later for the non- 
linear optical spectra, by Eq. (4.2), 



E e i*pm(g)g Q ei<p'm'(g)- (4.2) 



2Npp 
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Pa 
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The sum over the carbon_ppsitions can be evaluated using 
the symmetry properties^!!: 



5 pp ,V5D(ll\p -p) + 6 p ^ p ,D(22\-p p) 



(Ipm \ex a \ I'p'm') = 

5 p ^ p ,D(ll'\pp') £ (-^'C-m'^M, 
M—-1 

where the effective dipole moment is given by 
^(21 + 1)(2Z' + 1) 



: ( — 1) W E ^2NX-N' A aM' 
M" = -l 



(4.9) 



(4.3) 



with D(ll| - +) = D(ll| + -) = -0.70700ei? and 
D(22| + -) = D(22| - +) = 0.75585ei? - The dipole- 
allowed transitions from Txu to both Hq and Tig ex- 
citons are possible because of the difference of the Ceo 
symmetry group from the spherical symmetry. 



D(ll'\pp') = eR 



'pp'\ee)c x + Gio{ll'pp'\ee)e z , (4.4JB. Linear absorption: Optically allowed excitons 



with Gln{H'pp'\^) defined in Eq. ([T^) for the reference 
atom at position e, and the projection of the Cartesian 
components onto the L — 1 angular momentum states, 



Consider first the case of the Aq ground state. Its 
optical properties are governed by the time-dependent 
polarization field 



— ^ ax ~^ — &Ml) 

1 



P a (t) = 2n(0 P a (t) 



(4.5) 



The optical transition matrix elements from the Aq 
ground state to the electron-hole pair excited states in- 
troduced in Eq. (3.3) are 



(4.10) 



where n is the density of the buckyballs. The linear re- 
spon se i n the rotating- wave approximation is, following 
Eq. (4.6), given by 
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c Ap (l-)D(12\-p p) 
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0) = 5 L1 5 Pc ^ Ph D(12\p ePh )Al N . (4.6) 
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Selection rules limit excitations, to L = 1 excitons with 
odd parity, i.e. Txu statea!4E3. Out of the 60 singlet 
excitons there are only six such states, namely \1N + 
- > and \1N — h > with N = —1,0, 1. Because of the 
differ ent structures of the single-particle eigenvectors in 
Eq. (2.8) the oscillator strengths vary with the parity- 



(4.11) 



allowed pairs: 



£>(12| + -) = -0.48547ei? 
D(12| - +) = -0.55626ei? - 



(4.7) 



The characteristic dipole length is about half the distance 
from a carbon atom to the center of the cage. 

The coupling of two different excitonic states by an 
external electric field via the polarization operator is 



where u> is the frequency of light and a phenomenological 
lifetime-broadening parameter Tlp\ has been introduced 
for the electron-hole pair states. The diagonality and 
isotropy of the susceptibility tensor follows immediately 

from Em=-1 A aM A 0M = <W 

The resulting low-energy absorption spectrum of the 
7r-electron system is shown in Fig. 6 for various values 
of the effective intraatomic Coulomb interaction U. The 
ratio of the interatomic and intraatomic Coulomb inter- 
action has been fixed at the values e 2 /(eRoU) — 0.2377 
(left panel) and 0.0517 (right panel). Using an intrasite 
matrix element U = 17.3 eV the two values correspond 
to the cases of no screening of the interatomic interac- 
tion (e = 1) and of intermediate screening (e = 4.6). In 
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the spectral range considered the absorption in the tt- 
electron system exhibits two T\u exciton peaks at about 
w V + 0.012C7 - 0.666^- and w 1.139V 

+0.039J7 - 0.488-tJ- with different oscillator strengths. 

Comparison of the calculated absorption spectrum 
with experiment leads to a discussion of two salient 
points: (i) the comparison of the calculated and mea- 
sured peak positions and (ii) the relative intensities of the 
two absorption peaks. The two calculated exciton peaks 
should be compared to the measured absorption peaks at 
3.81 and,4.90 eVEJ, or 3.65 and 4.72 eVEl, or 3.78 and 
4.84 eVclL The weak structure observed at lower energies, 
in particular the weak peak at 2.73 eV, for solid CgocHl has 
been ascribed to dipole-forbidden transitions, which be- 
come partially allowed due to lattice fluctuations, inter- 
face effects and/or internal electric fields. The same holds 
for the small structure occurring in the energy region of 2 
eV in spectra of Ceo isolated in a noble gawpatrixEE TV 
intense absorption band close to 5.96 eVc2l or 5.87 e'\ 
only possesses a partial n character and is also beyond 
the scope of our study. The identification of the lower 
and upper T\u excitons (L = 1, P = — ) with the two 
absorption peaks under consideration restricts the range 
for the values of the parameters U and e 2 /(eRo). If we 
take the strong onsite interaction value U = 17.3 eV, the 
dielectric constant varies drastically between e w 2 (for 
the hopping parameter V = 4.62 eV) and e w 10 (for 
V = 3.83 eV). With a weaker U « 10 eV, the range of 
the dielectric constant is reduced to 2.5 < e < 4.8. 

The theoretical result that the absorption peak at 
lower energy possesses a smaller oscillator strength than 
the peak at higher energy is in agreement with the exper- 
imental observation. Within the single-particle approxi- 
mation, i.e., for U = and e 2 /(ei?o) = 0, the oscillator 
strengths of the two transitions h u — * t\ g and h 



in 



between molecule levels are similar, from Eq. (4.7). The 
reduction of the ratio of the oscillator strength of the 
lower-energy peak to that of the higher-energy peak may 
be understood in terms of the Coulomb coupling of the 
two Tm excitons, as is evident fr om t he presence of the 
eigenvectors CA P (1-) from Eq. (fjj|) in the oscillator 
strength of Eq. ( 4.1l|) . Since the dipole matrix elements 
Z?(12| H — ) and D(12| — h) have the same sign and nearly 
equal magnitudes, the oscillator strengths of the two T\u 
excitations are governed by the relative sign of the co- 
efficients ca p (1— ), and so by the sign of the coupling 
term V-\ (1— ) in Eq. ( 3.10] ). The domination of the ex- 
change term leads to V] (1— ) < and, hence, the oscil- 
lator strength of the low-energy exciton at Hlu = E\ 

is reduced compared to the high-energy absorption at 

Tvjj = Ei |_. This shows not only that the relative 

strengths of the two exciton peaks are influenced by the 
Coulomb interaction but also that the intersite interac- 
tion is indispensable. 

Let us consider briefly the consequences on the lin- 
ear opt ical spectra of the fascinating possibility raised 
in Sec. Ill that the excitonic state T 2 g is the ground 



state rather than the closed-shell Aq state. In the en- 
ergy region of the two allowed single-particle transitions 
h u — > tig and h g — > ti u , two peaks can arise from the 
transitions between the lowest T 2 g state to the two Hu 
and two Gjj states (see Fig. 4) satisfying both the angu- 
lar momentum and parity selection rules. Their oscilla- 
tor strengths are governed by the moments -D(ll| H — ) 

and £>(22| h). The transition energies may also be 

accounted for by the parameter values for V, U, and 
e 2 /(ei?o) already discussed in connection with Fig 5. 
However, the single-particle energies from the new exci- 
tonic ground state have to be recalculated and compared 
with the measurements by photoemission and inverse 
photoemission. Since in the exciton ground state there is 
already one electron in a LUMO and a hole in a HOMO 
state, the one-particle excited state has strong correlation 
effects. The interaction effects can split the degeneracies 
of the HOMOs and the LUMOs, leading to more than 
two electron and two hole single-particle states. The 
group-theoretical indentities t\ p x Tig — 9p + hp an d 
h p xT 2 G = tip+hp+g P +hp (p = g,u) give us a rough idea 
of the single-particle multiplets. The symmetry breaking 
which will remove the degeneracy of the Tig states would 
reduce the number of one-particle excited states. 



C. Electro-optic Kerr effect: Electric-field-induced 
forbidden excitons 

In the theoretical linear optical spectra, the parity se- 
lection rule excludes the same-parity transitions, such as 
those lowest in energy, h u — > t\ u . They can, however, be 
induced by the application of a static electric field which 
mixes the even-parity excitons with the odd-parity ex- 
citons. The linear response t o the external laser field of 
the electric polarization, Eq. ( 4.10| ), to second order in an 
applied static electric field F (the Kerr effect), is closely 
relat ed to the third-order susceptibility, to be deri ved i n 
Sec. [VP . By a similar derivation leading to Eq. ( 4.15 ), 
the static electric field effect on the linear optical response 
is given by 



F 2 



E E 

£=1A', A"=± 



(4.12) 



St, A' S LA'' 
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with the prefactor depending on the direction of the ap- 
plied static field denoted by its unit vector F relative to 
the light polarization direction: 



Pai3(F)^—(3S al3 + F a F ), 



(4.13) 



and the oscillator strength Sla given in Eq. (4.17). 

The selection rules AL — 0, ±1 confine the contribu- 
tions to the Kerr effect from the Hg and Tig excitons. 
Because of the random orientations of the Cqq molecules 
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in solutions or in the face-centered cubic thin films 
the dependence of the field direction of the prefactor is 
in practice averaged out to P a p{F) = 5 a p. Only in the 
case of the low-temperature simple cubic crystals of un- 
doped Cqq is there a chance of experimentally partially 
probing the polarization dependence. The four molecules 
per unit cell are rotated by an angle <f> around a space 
diagonal axis [111], [111], [111], or 13X11 The angle of 
rotation is found to be (j) = 22 - 26 El~y. 

The spectrum of the electrooptic Kerr effect is pre- 
sented in Fig. 7 as a function of the reduced pho- 
ton energy for different model parameters of the effec- 
tive Coulomb interaction, U and e 2 /ei?o- The electric 
field strength has been fixed at eFR /{2V) = 0.001 
to compare with the allowed transitions in zero field 
in Fig. 6, corresponding to 2.6xl0 5 V/cm for a hop- 
ping parameter of V — 4.626 eV. The positions of 
the field-induced excitons are approximately given by 
£ 2 +- = 0.757V + 0.019E/ - 0.760e 2 /ei? for the H G ex- 
citon or E 1+ _ = 0.757V + 0.012117 - 0.986e 2 /ei? for 
the T\g exciton. Their energy difference is determined 

by the interatomic Coulomb interaction, i?2H — E\-\ w 

0.226e 2 /ei? ~ 0.93 eV/e. For intermediate screening, it 
is smaller than the line broadening used in Fig. 7, explain- 
ing the appearance of only one pronounced field-induced 
peak in the absorption spectrum below the energy of the 
lower allowed T\jj exciton (in Fig. 4). Only for very weak 
screening does a weak second Tig-related peak appear at 
an energy below the more intense -He-related peak (see 
the solid and dashed curves in Fig. 7a). 

These general findings seem to be in agreement with 
the results inferred from different nonlinear optical ex- 
periments, such as two-photon absorption^ and degen- 
erate four wave mixingo. A two-photon resonance at 
2.73 eV or 2.67 eV is observed and identified with an Hq 
exciton. Such energy values arc within the range of the 
theoretical values given by the parameters under discus- 
sion. The prediction of a field-induced Hq exciton line 
alsp-|explains the weak absorption structure around 2.73 
eVc2l. The theoretical energy ordering T\q < Hg < T\u 
clearly visible in Figs. 4 and 5 agrees with results, 
the second-harmonic generation (SHG) on a surfacecZH 
where a resonance near 1.8-1.9 eV is identified with a 
Tig exciton. However, the theoretical prediction of the 
energy splitting between Tig and Hq is smaller than the 
experimental value of 0.8-0.9 eV if we choose our model 
parameters, such as U = 17.3 eV, e — 4.6, to obtain 
similar exciton energies as the quantum-chemical calcu- 
lations. 



D. Third- harmonic generation 

In this section, we consider another nonlinear opti- 
cal experiment, the third-harmonic generation (THG). 
A three-photon resonance occurs when three times the 
fundamental photon energy is equal to the lowest Aq — > 



T\u one-photon dipole-allowed transition (in the single- 
particle picture: h u — > t\ g and h g — > From the gen- 
eral expression for. the third-order susceptibility given by 
Armstrong et alo, we take into account only triply reso- 
nant contributions to the susceptibility describing THG. 
Neglecting biexciton effe ctsEj and using the dipole matrix 
elements (fh|) and (4.9) we find 



2 +L +1 

= i^E E E E 
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Ps 
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(4.14) 



Here, n is the density [cf. Eq. (4.1C )] and P/3 7 5 denotes 
the sum over all permutations of /3, 7, and 5 ensuring 
that the fourth-rank tensor third-order susceptibility is 
independent of the ordering of those three indices, i.e., 
the Cartesian components of the three fields creating the 
third-order harmonic. 

Selection rules dictate that the three-photon resonance 
takes the system from the Ag ground state to the L = 1 
excitons with odd parity (Tit/), and from T\u to L = 2 
(Hq) and L = 1 (Tig) excitons with even parity, which 
may give rise to the doubly resonant terms. We restrict 
ourselves to the double-resonance terms with the lower- 
energy (A = — ) excitons. The higher- lying excitons are 
energetically well separated from the frequency region of 
interest (cf. Fig. 4). Consequently the \LN H — ) states 
may be approximately replaced by the uncoupled \LNpp) 
(p = +, L = 1,2) states, corresponding to the single- 
particle transitions h u — > t\ u . The third-order suscepti- 
bility becomes 



L = l A',A"=± 



°LA' 

(£i_A' - 3ftw - £Ti_a') 
Sla 1 ' 



(4.15) 



(El+- - - ir L+ _) (-Ei_a» - hu; - iTi-A'i) 
with the polarization-dependent prefactor 

1 +2 

G 0/J7 a = -P^, Kfi{M)F 5 ^(M), (4.16) 
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+ 1 

F af3 (M)= £ (-l) M 'Ai MI C^';_ M ,A^ MII , 
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and the oscillator strength 



P =± 
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Eq. (4.15) may be interpreted as the third-order suscep- 
tibility of a five-level system: ground state Aq, even- 
parity excited states Hq or Tig, and the two odd-parity 
excited states T\u , in contrast to the third-order suscep- 
tibility of a three-level system used to fit the experimen- 
tal measurementsLal£3. Moreover, the symmetry of the 
states participating in the THG process are well defined 
here. On the other hand, only the most dominant res- 
onant terms are considered here, whereas the formulas 
used to fit the data also include less important resonant 
and nonresonant terms. 

The spectral behavior of the magnitude of the third- 
order susceptibility is plotted in Fig. 8 for the same set 
of model parameters as in Fig. 7. Only the lower T\u ex- 
citons appear. Despite the presence of many exciton lev- 
els, in the spectral region considered only a pronounced 
double-peak structure or one broad peak appears. In the 
case of stronger screening of the intersite Coulomb inter- 
action (the right panel of Fig. 8), the high-energy peak 
corresponds to the two-photon resonance with the L = 2, 
P = + (Hg) exciton. It is enhanced by the resonance 
with the Tig excitons (L — 1, P — +). The three-photon 
resonance at the frequency of the fundamental light wave 
occurs at slightly lower energies defined by the electric- 
dipole-allowed L = 1, P = — {T\u) exciton. On the 
other hand, in the case of the weaker screening (the left 
panel), for the curves going from the right to the left in 
order of increasing onsite Coulomb interaction U , there 
is an interchange of the energy order of the two reso- 
nances hu) = \E\ and Tiuj = \E 2 ^ , \E\ + ~. The 

near coincidence of the triple and double resonances for 
the intermediate values of U creates the appearance of 
a strong single peak in Fig. 8a. For larger values of the 
onsite Coulomb interaction, the Hg and Tig resonances 
split again and a weak Tic-related peak occurs at the 
low-energy tail of the THG structure. 

Our calculated spectrum can be used to infer the sym- 
metry of two-particle elementary excitatipas-observed ex- 
perimentally. In the THG experimentOO two peaks 
were observed at 1.3 /im and 1.06 fim. By taking a 
combination of relatively large parameters V , U, and 
e 2 /(ei?o), we can interpret the lower two-photon reso- 
nance to yield the measured energy of the Hg exciton 
at 1.9 eV and the higher triple-photon resonance to yield 
the energy of the T\jj exciton at 3.5 eV. A previous inter- 
pretation of the low-energy peak in the THG spectrum 
as a two-photon resonance with the one-photon forbidden 
Tig levelEil, even though in agreement with the low value 
of the resonance energy, could not explain the absence 
of the Hg exciton in the spectrum. We have demon- 
strated here that, with the help of a careful symmetry 
analysis, nonlinear optical spectroscopy can clarify the 



complicated electronic structure of the Cqq molecule, es- 
pecially its electron-hole pair excitations. 

We note briefly that, in the case of the excitonic ground 
state T2G, the interpretation of the THG would be dif- 
ferent. According to the AL = ±1 and AL = se- 
lection rules, a three-photon resonance excites the Hjj 

or Gu exciton at Tiuj = (E L - £? 3+ _)/3, (L = 2,3). 

Many transitions may contribute to the two-photon res- 
onance. Among them are Hu — > Tig, Hg, T 2 g, and Gg, 
or Gu — ► Hg, Tig, and G G . 



V. DISCUSSIONS 

In this paper, we address the issue of the strong elec- 
tron correlation in a Cqq molecule by an analytical model 
making maximum use of symmetry. The basic compo- 
nent is a simple quasi-particle model based on a tight- 
binding scheme of 7r— orbitals of the HOMO and LUMO 
states, with parameters deduced from photoemission and 
inverse photoemission experiments. The assumption of 
the closed-shell ground state is tested against the energies 
of the electron-hole pair excited states. The interaction 
between electrons and holes includes the Coulomb inter- 
action of the 7T— orbitals on the same carbon atom and 
the long-range interaction between different carbon sites 
with a constant dielectric screening. The onsite electron- 
hole interaction is dominated by the exchange term, and 
is, therefore, repulsive while the electron-hole interaction 
on different carbon sites is attractive. The balance be- 
tween these two terms depends, in our model, on the 
dielectric screening. 

Our theory of the linear optical absorption spectra 
based on the closed-shell ground state agrees with ex- 
periment for a reasonably strong dielectric screening of 
the long-range interaction. We have established the rela- 
tion between the energies or relative oscillator strengths 
of the dipole-active excitons and the Coulomb interac- 
tion terms connecting these excitons of the same sym- 
metry. This coupling effect explains the observed strong 
T\u exciton peaks in the linear optical absorption spec- 
tra. For the third-harmonic generation, our theory iden- 
tifies the doublet excitation at the fundamental wave- 
length A w 1.3 /im as the forbidden Hg exciton and not 
as the Tig exciton suggested by experiment, whereas the 
triple resonance is clearly related to the T\u excitation. 
This identification is supported by the argument invok- 
ing the approximate angular momentum selection rule in 
the nearly centrosymmetric molecule (see Sec. IV). The 
same transition to the Hg pairs is found to play a role 
if an external static electric field is applied to the Cgo 
molecules. We predict a pronounced optical Kerr effect 
with a photon resonance at the Hq energy. 

We wish to note an interesting scenario in which the 
closed-shell Aq state is unstable against an electron-hole 
pair excitation of symmetry T 2 g with a reasonably weak 
dielectric screening. The consequences of having a T 2 q 
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exciton state as the ground state are extraordinary. The 
optically dipolc- allowed states are subsets of the Hjj. The 
most interesting possibility is that a quasi-particle state 
now consists of an electron (or hole) plus an electron-hole 
pair, leading to strong correlation effects. Symmetry con- 
sideration would lead to more photocmission lines than 
the closed-shell ground state scenario. An experimental 
test of the electro-optic effect, which can discriminate the 
ground state symmetry, is described in Sec. IV. C. 

A less extreme possibility is the small excitation en- 
ergy of the T 2 q exciton. The quasiparticle dynamics can 
be affected by the easy Coulomb excitations of such exci- 
tons. In particular, this could provide a source of effective 
quasiparticle interaction and, thus, possibly a source for 
superconductivity. 

While our model study of the excitons has indicated 
that the most likely scenario of the ground state in Ceo 
is the closed-shell Aq, it points out the possibility of 
constructing other molecular solids which lowers the ex- 
citonic states. The requirement is the strength of the 
inter-site electron-hole attraction over the onsite repul- 
sion, which might be achieved by increasing the number 
of atom sites in a molecule. Our model approach may 
also be used to study excitons in quantum dots when the 
effective-mass approximation fails. We leave the investi- 
gation of the above conjectures to future work. 
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FIG.l: The structure of the Cqo molecule. Each dot 
denotes the position of a carbon atom. The connecting 
solid lines indicate the bonds. The intact icosahedron 
with the vertices at N, S, A i7 and A[ (i = 1, . . . , 5) is 
indicated by the dashed lines. 

FIG. 2: Schematic four-level diagram of HOMO and 
LUMO states for characterization of the lowest electronic 
excitations in the Ceo molecule. 

FIG. 3: Feynman diagrams for the electron (upward ar- 
rows) - hole (downward arrows) interaction (dashed line): 
(a) the direct attraction, (b) the exchange counterpart. 



FIG. 4: The excitation energies of Frenkel excitons be- 
longing to the single-particle pairs h u — ► t\ u , h u — > t\ g , 
h g — > t\ u , and h g — > t\ g versus the effective interatomic 
Coulomb interaction strength e 2 / (ei?o)- The intraatomic 
Coulomb interaction is fixed at two values (a) U = and 
(b) U = AV. Tip: solid line, H P : dashed line. Both 
parities P = +, — arc considered. The high-energy (low- 
energy) solutions of the coupled pairs A = + (A = — ) 
are plotted as thick (thin) lines. All energies are given in 
units of the hopping parameter V. 

FIG. 5: The lowest even-parity exciton energies versus 
the interatomic Coulomb interaction. The zero line is 
given by the closed-shell Hartree-Fock ground state, (a) 
U = 0, (b) U = AV. 

FIG. 6: Absorption spectra of 6*60 versus the photon 
energy in units of V for different parameters of the ef- 
fective Coulomb interaction and damping parameters T. 
Solid line: U = 3V, dashed line: U = 2V, dotted line: 
U = V, dash-dotted line: U = 0. The ratio of the 
intersite and intrasite Coulomb interaction is fixed to 
e 2 /(eR U) = 0.2377 (left panel) and 0.0517 (right panel). 

FIG. 7: Absorption spectra of Cqo molecules in a static 
electric field versus the photon energy near the forbidden 
transition h u — > t\ u for different parameters of the ef- 
fective Coulomb interaction. Solid line: U = 3V, dashed 
line: U = 2V, dotted line: U = V, dash-dotted line: U = 
0. (a) e 2 /{eR U) = 0.2377, (b) e 2 /{eR U) = 0.0517. The 
damping parameter is chosen as T = 0.05V. Each spec- 
trum has to be multiplied with the prefactor (eFR /2V) 2 
to compare with the strength in Fig. 6. 

FIG. 8: Spectral variation of the THG susceptibility 
in the region of the lower 3-photon Tm and 2-photon 
Hq resonances for the same parameters of the effective 
Coulomb interaction as in Fig. 7. The damping parame- 
ter is chosen as T — 0.031^ for both allowed and forbidden 
excitons. 
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Table I. The intraatomic and interatomic Coulomb interaction matrix elements ( p. Ill , p. 12j , p. 13| ) . 
(I.R. — Irreducible representation) 



I.R. 




Quantum numbers 




F ppl (LNP) 






H PP >(LNP) 






X ppl {LNP) 


j 




L 


JV 


P 


+ + 




+- 


++ 




+- 


++ 




+- 


Tig 


1 


0,±1 


1 


0.50000 


0.76690 


0.61923 


0.87251 


0.89689 


0.22713 


-0.02111 


-0.03239 


-0.02615 


H G 


2 


0,±1,±2 


1 


1.50000 


1.23310 


0.79877 


0.85182 


0.82745 


0.05389 


0.03329 


0.04047 


0.05790 


Tig 


3 


0,±3 


1 


1.33333 


1.99975 


-1.63289 


0.85527 


0.87308 


-0.09762 


-0.08122 


-0.12181 


0.09947 


Gg 


3 


±1,±2 


1 


0.50000 


0.13364 


-0.23822 


0.87251 


0.87134 


-0.16449 


-0.00908 


-0.00303 


0.00416 


Tiu 


1 


0,±1 


-1 


1.54204 


1.50000 


1.50777 


0.85363 


0.84560 


0.19439 


0.16148 


0.11124 


0.13512 


Hu 


2 


0,±1,±2 


-1 


0.45796 


0.50000 


-0.08977 


0.87070 


0.87873 


0.08662 


-0.01309 


-0.02673 


0.00445 


T2U 


3 


0,±3 


-1 


0.63864 


1.00000 


0.21489 


0.86786 


0.91547 


-0.18339 


0.02377 


-0.01409 


0.01251 


Gu 


3 


±1,±2 


-1 


1.54204 


1.25000 


-1.17978 


0.85363 


0.81391 


-0.11652 


-0.06371 


-0.07559 


0.07082 
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E 1 _(l1m->) 

•H- -H- -H- "H- 4t E2-(l2m->) 
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